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SUMMARY 


W 

1 The directional and frequency diffusion of a plane monochromatic 

2 sound wave in statistically homogeneous, isotropic, and stationary tur- 

7 bulence is analyzed theoretically. The treatment is based on the diffu- 

sion equation for the energy density of sound waves, using the scattering 
cross section derived by Kraichnan for the type of turbulence assumed 
here . 


A form for the frequency -wave number spectrum of the turbulence is 
adopted which contains the pertinent parameters of the flow and is 
adapted to ease of calculation. A new approach to the evaluation of the 
characteristic period of the flow is suggested. This spectrum is then 
related to the scattering cross section. 

Finally, a diffusion equation is derived as a small-angle scattering 
approximation to the rigorous transport equation. The rate of spread of 
the incident wave in frequency and direction is calculated, as well as 
the power spectrum and autocorrelation for the wave. 


INTRODUCTION 


The fluctuation of acoustic signals in media of randomly varying 
index of refraction has been considered by several investigators from 
nearly as many points of view. On the basis of geometrical or ray acous- 
tics, Bergmann (ref. l) and Obukhov (ref. 2 ) have evaluated the fluctua- 
tions in phase and amplitude of a spherical acoustic wave emanating from 
a point source. Skudrzyk (ref. j) has also estimated these quantities 
and obtained the probability distribution of the signal from a geometri- 
cal scattering point of view and has tabulated his results. Mintzer 
(refs. 4, 5 , and 6 ), using single and double scattering formulations, 
also calculated similar quantities. One of the beauties of Obukhov's 
approach is that he was able to obtain both Mintzer' s and Bergmann 's 
results as limiting forms of a single derivation. 
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The present work proposes to solve somewrat the same problem from 
a multiple scattering point of view by formulating the transport equa- 
tion for the acoustic energy density spectrum in turbulent flow. The 
diffusion takes place both in frequency and direction, since the flow 
varies in time as well as space. The directicnal diffusion of a sound 
wave has been previously considered by Lighthill (ref. 7)* T ^ho proceeded 
from a diffusion equation which appears to have been inspired by random 
walk considerations . 

The present investigation was conducted at the Acoustics Laboratory 
of the Massachusetts Institute of Technology under the sponsorship and y 

with the financial assistance of the National Advisory Committee for ^ 

Aeronautics. 2 
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SYMBOLS 


A,B parameters for hypothetical spectrum 

b rms frequency variation 

c wave in undisturbed medium 

dk volume element in k space 

die volume element in x space 

a(s>) element of solid angle 

E wave number and spatial density of acoustic energy 

£ acoustic energy 

J 2 acoustic intensity 

I wave number and spatial density of acoustic intensity 

I intensity of unscattered sound wave 

J Q (x) Bessel's function of zero order 

K scattering vector and turbulent wave number 

k wave vector 
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wave number where W(K) has its maximum 
correlation length 

Mach number of turbulent fluctuations 
summation index 

unit vectors in direction of k and k T , respectively 

Legendre polynomial 

power spectrum of scattered wave 

total scattering cross section per unit volume 

frequency variable 

time 

space -time transform of u 
turbulent velocity field 
volume of scattering region 

isotropic spectrum for time -independent turbulence 
isotropic spectrum for time -dependent turbulence 
time transform of acoustic condensation 
position vector 

distance over which sound has propagated in atmosphere 
inverse of k Q 

scale factor in W(K) 

gamma function 

normalizing constant 

relative frequency shift from tx> Q 

Dirac delta function 

dissipation rate for turbulent energy 
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v 

parameter in hypothetical spectrum 


n 

distance variable 


9 

total angle of deviation from original direction 


K 

transform variable 


A 

characteristic frequency of turbulence 


V 

kinematic viscosity 

W 

1 

cos 0 

1 

2 

P 

relative frequency shift from cu 

7 

Po 

ambient density of atmosphere 


p 

branch points of integrand 


1 



2 



T 

delay time 

- 

a 

differential scattering cross section 


a 

vector from x r to x 


0 

scattering angle 


<fc. . 
10 

spectral function for homogeneous tirbulence 



azimuthal angle 


'K'O 

autocorrelation function for acoustic signal 


ft 

frequency of turbulence 


O) 

frequency associated with wave number k 



frequency of unperturbed sound wave 


Superscripts : 


( )■ 

pertaining to scattered quantity 


C) 

vector 
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Subscripts : 
' >av 

<•>„ 


P 


ensemble average 

averaged over differential cross 
propagation 


section 


o 


TRANSPORT EQUATION 

The spectral energy density E(x,k) of a sound field is defined in 
such a manner that in the volume dx and in the range of wave numbers 
dk, (i.e., the wave number k falls in the range k, k + dk) the acous- 
tic energy is 

d t = E(x,k)dx dk (l) 

A beam of intensity d J and the wave number in dk has energy scattered 
into beams with wave numbers in dk’ at a rate given by | du?| a(k/k 1 ) dk T . 
In the plane waves being considered, the intensity is related to the 
energy density by 

dj = c — E(x,k) dx dk = l(x,k)dx dk (2) 

k 

where l(x,k) is an "intensity spectral density" and c is the wave 
speed in the undisturbed medium. It is assumed that the rate of change 
of energy in time is totally accounted for by scattering and propagation, 
and such effects as viscous, thermal, and relaxation loss mechanisms are 
neglected. 

The energy loss rate due to propagation alone is given by 


QE(x,k) 


dt 


= -div I 


( 3 ) 


In the volume dx, the conversion of k waves to k' waves results in 
dE(x,kjl 


dt 


Jk-> k 


_ = - f | l(x,k)| a(k/k’ )dk' 
k r J it t 


( 4 ) 



and the conversion of k 


waves to 


k waves gives 


<3E(x,K) | 

- ' k ' k 



|l(x,k') 0(k'/k)dk' 


( 5 ) 


The integration in equation (4) may be perfoimed immediately if the total 
cross section is introduced 



a(k/k' )dk ! 


( 6 ) 


If the acoustic source has come into statistical equilibrium with the 
field, then the total rate of change of energy vanishes. Thus, adding 
equations ( 3 ), (4), and ( 5 ) and setting the result equal to zero gives 
the transport equation 


div I = -cE(x,k) + c / dk'a(k '/k)E(x,k' ) (7) 

J k* 

Now it is necessary to evaluate o(k/k') for turbulent scattering. 

Scattering Cross Section a(k/k') 

The expression for the time -dependent cross section used here is 
essentially that derived by Kraichnan (ref. £) . Let the directions of 
k and k' be defined by the unit vectors n and n', respectively. 

The turbulent velocity field is described by a Fourier transform, 

U(K,H) = dxj" dte- i(< ' x - nt) u(x,t) ( 8 ) 

Referring to appendix A and reference 8 , the reader will see that the 
appropriate form for a is 

o(iVk') = '^f(n ■ fl’) 2 (|l„| 2 ) aT ( 9) 

where U n is U • n and o> is the frequency of the incident sound. The 

ensemble averaged factor in equation ( 9 ) must be related to the one- 
dimensional spatial spectrum for isotropic turbulence. This spectrum, 
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as described by Batchelor (ref. 9), is denoted by W(k) . The kinetic 
energy and the rate of energy dissipation per unit mass are (ref. 9) 

p 00 

\ (u 2 ) = / W(K)dK (10) 

^ Jo 

and 



For time -dependent turbulence, the spectrum w(K,n) is introduced which 
must satisfy 


W 



w(K,fi)dn 


( 12 ) 


The mean flow is assumed to vanish and the autocorrelation at a point is 
assumed of the form 


<t 1 (t)u 1 (t + = ^u 2 (t)u 2 (t + t)^ 

+ t)^> 


- \ J 2 \ e-AUI 

3 \ /av 


(13) 


with a frequency spectrum 


3 ^ jt(a 2 + fi 2 ) 


If it is assumed that the frequency and wave number spectra may be sepa- 
rated, the equation can immediately be written 


v(K,n) = w(K)— — £ v 

*(a 2 + ft 2 ) 


(14) 


where equations (l2) and ( 13 ) were used. A convenient form of W to 
choose for the integrations to be performed is 


W(K) = apK^e -1 ^ 


(15) 
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where l/a = k Q is the wave number for the maximum of the spectrum; 
that is, it is the wave number of the energy -bearing eddies. This spec- 
trum is chosen since it has the known behavior at small wave numbers 

(ref. 9) with the energy -bearing maximum and allows analytic evaluation 
of the integrals to follow. Using equation ( 10 ), 

3 > ^<u 2 > (16) 

40 

which ties the parameters a and |3 to important quantities in the 
description of the turbulent field. 

It is shown in appendix B that ^U n j^ may be related to the spec- 
trum equation (l 4 ) by the relation 

(l U n^ a )] 2 ) av '^| Zl “ (K ' I1 > 


where 0 


is the angle between incident and scattered wave vectors. 


The characteristic frequency of the turbulence by equations (13) 
and (l 4 ) is A. This may be approximately e/aluated from parameters of 
the flow by setting 


A 



(18) 


where u = \j \u / and L is the correlatioi 
(mean eddy size). Lighthill (ref. 7) gives 1 
which is used in this investigation. It is 


length of the turbulence 
formula for this length 


L = 3* f°° W(1 0 ^ . 3*a 

2 < u 2 >-'0 K " 4 

for the spectrum equation (15). 


(19) 


There is another possible method for thi evaluation of A. In an 
equilibrium turbulent field there is a balance between the^ randomly 
created vorticity and its dissipation as viscous loss. This is reminis- 
cent of the equilibrium state of a harmonic oscillator (with damping) 
when excited by random noise. For the oscillator, one obtains a correla- 

^ 1 | ^ i 

tion function ~e where A is given by the ratio of the energy 

lost per second to twice the average energy stored. Making this analo- 
gous to the turbulence problem, 


:* h oj r- 



2T 
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In terms of the spectrum equation (15), this is 


( 20 ) 


A = 


15 v 
8a 2 


( 21 ) 


which is lower by a factor of 10 or so than experimental results which 
support the value given by equation (l8). It is believed that the fail- 
ure of equation (21) is due primarily to equation (15) and not to the 
idea expressed by equation (20) . The assumed spectrum cuts off too 
sharply for large wave numbers and hence gives too low a value for e. 

Any spectrum which gives the range predicted by equilibrium 

theory, however, appears to be too complex to integrate in expressions 
which arise later in the paper. Hence, equation (15) will be used recog- 
nizing its limitations at large wave numbers. 


where 


k 


A possible spectrum is W ~ K^l + AK ^/3l^ + BK 

p < 0 but arbitrary. ‘This spectrum contains the necessary K , , 

and K”^ ranges, but integrations involving this factor are seen to be 
rather formidable. 


K 


■5/5 


DERIVATION OF DIFFUSION EQUATION 


Referring now to equations (9), (l4), (15), and (17) and to 
reference 8, 

2 _ -1 -haK 

a( k/k 1 ) = 03 a Y? — (n * n' ) 2 cos 2 (0/2) (22) 

4c 5 A 2 + Q. 2 

where K = |k| = |k - k'| and Q = j cu - o>'| . The frequency shift is 
and the magnitude of the scattering vector K. It is seen that 


o(k/k') » a(k'/k) 


(23) 


if the change in frequency upon scattering is not too great. 

An attempt was made to simplify equation (7) by some appropriate 
assumptions. If the incident acoustic wavelength is smaller than the 
eddy size, then the scattering will be strongly forward. Also, the 
period of the sound wave is usually much less than the period of the 
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turbulence (oj » Q). Thus, a(k/k') is nonzei o for only small values 
of k - k‘, and E(x,k') can be expanded arour d (x,k) in a Taylor's 
series . 

Since the diffusion of a plane wave of frequency o> 0 with a wave 
vector parallel to the z-axis is to be considered, it is expected that 

CD - CD r 

the function E will depend only on z, B, ard A = — y where 8 

is the polar angle from the z-axis. If \ = cc s 8 , then 




where /f(| , ,A , )\ = j f (| 1 ,A ! )a(k/k 1 )dk 1 . Putting equation (24) into 

x 7 a uk' 

equation ( 7 ) > gives 



W 

1 

2 

7 


This is the diffusion equation which will be used. With a knowledge 
of o y it remains to calculate the unknown coefficients. 

If p = A f - A, using the law of cosines 

K = k 2 p 2 + 4k 2 sin 2 (0/2) - 4k 2 p siri 2 (0/2)J ( 26 ) 

Now examine the relative magnitudes of these terms. Consider the ratio 
kp/2k sin(0/2) . Now kp « A/c = AM/u, where A is a characteristic fre- 
quency of the turbulence and M is the Mach number of the turbulent 
fluctuations. The factor 2k sin(0/2) « 1 fl y wiere l is a character- 
istic length for the flow. Then 

kp _ AM 
2 k sin(0/2) u7* 


M « L 


(27) 



^ H OJ ^ 
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substituting for equation (l 8 ) . Also, p « A/cd 0 « 1, This means 
equation ( 26 ) can be expanded as 


K * 2kx(l - A + A P 2 
\ 2 8 x 2 / 

and x = sin(0/2) . Using equation (28) in (22) gives 
a(k/k') ~ ^\: J+ x 2 (l - x 2 )(l - 2 x 2 ) 2 (A 2 + a> 2 p 2 ) -1 exp 


( 28 ) 


-8akx jl 


* + ^ 


(29) 


COEFFICIENTS IN DIFFUSION EQUATION 


To calculate the "averaged" coefficients in equation (25) consid- 
ering ' -i) 

(I ' - S / 0 = J dk’a(k/k' ) ( £ ' - |) 

The evaluation of this integral using equation (29) is carried out in 
appendix C. Since | is independent of ft, it is possible to calculate 
this average with the simplified spectrum equation ( 15 ). The result is 
shown to be 

< s ' - 5 >*"(gr ) 5 (50) 

where = (u^/c^. 

The calculation of 1 - £)^) 0 is performed in appendix D, and again 

the simpler spectrum equation ( 15 ) can be utilized. This simplification 
occurs because w(K,ft) is chosen to be a product of 2 factors, one inde- 
pendent of K and the other, of ft. (it is probable that the actual 
isotropic spectrum which would be determined from the equations of motion 
would not factor in this manner, since it is seen from equation (l8) 
that the periods and eddy sizes are related.) The result here is 




(31) 
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The evaluation of (p) is not so straightforward as in the pre- 
vious cases. However , there is no physical he. sis for assuming that 
there will he any change of frequency of the vave. In appendix E, it is 
shown that although the expression (22) does predict some slight average 
frequency shift, this shift will be very small, and can he neglected if 
compared to its rms value. Accordingly, let 


(p) 0 = 0 (32) 

which also premits the assumption that ^p(§ ' - (■ )/ a = 0* 

The evaluation of j can carr f ec ^ forth without difficulty, 

and in appendix F it is found to he 


/p2 ) - a „2 

\ / a p^2 c 


(33) 


Notice that according to this result the rms frequency deviation should 
he proportional to A 1 / 2 instead of A as oi.e might suppose. 


W 

1 

2 
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SOLUTION OF DIFFUSION EQUATION 


By putting the calculated coefficients into the diffusion equa- 
tion (25), after making the proper definitions, one obtains 


dE = dE + l/i 2) Ae ^ 1 

dtl di 2 V '^ 2 2 dr 2 


(34) 


where r, = L, r = A and 7 — = H A L OA I 
32a 7 2 L5,j2c 

This equation must he solved with the condition that there Is a mono- 
chromatic plane wave at z = 0. The energy density at z = 0 will be 
of the form E = B&(l - |)6 (r), where B is 1.0 he determined. Then, 
from equation (l), the energy density becomes 


l (z = 0 ) = f E(0,k)dk = 2«Bk 0 5 r d|5(:. - i) f d(ry)5(r) 

J J 0 J -00 

or, I 0 C_1 = rtk Q ^7B, where I 0 is the intensity at 2=0. Then, 
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I n 

E(0,k) = s(r)B(l - 4) ( 35 ) 

nk^yc 

Equation ( 34 ) can be solved rather easily if the distances z are such 
that the initial wave has spread only slightly in direction. In this 
region then, let the term on the left-hand side of equation (3I+) reduce 
to c 3 E/drj. The solution in then immediately 


W 

1 

2 
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E( 5 ,K) = ^ J dKe“ 1Kr P m (4)A In (K)e 

m=0 


- T , ( v + ^ 111+1 i ) 


(36) 


The b 


oundary condition equation (35) defines A^( k) to be 


= p 1 + 1) 
4 rt 


(37) 


Finally then. 


/ CO 

E(-l,S,r) = — iS— H /2 e- A / 2, > y(2m + 

4 it^k 0 3 rc / 


(38) 


Z_j 

m=0 


which reduces to Lighthill’s result (ref. 7) for time -independent scat- 
tering. The separation of the frequency and angular diffusion processes 
into factors in equation (38) is, of course, a result of the form of the 

turbulent spectrum and the vanishing of the term /p(£ 1 - §)\ . 

/ a 

This equation may be made more transparent by considering the dif- 
fusion to have progressed only slightly. Then the angular deviations 
are small and it is possible to substitute 

— (l P 2 -)— d 2 1 
<3| ' " ^d4 ^ d0 2 + 0 §0 



for which the general solution is 


p 00 

E(x,k) = / pdpJ o (p0)e 

Jo 


\ P 2 4 


j: 


1 2 

litr - t| 

dKe A(k,P) 


( 40 ) 



where A(k,3) Is evaluated by requiring that E(x,k) correspond to a 
plane monochromatic wave at z = 0 of intensity I 0 . The result of 

this gives 

E(x,k) = I (3dpJ 0 (pe)e ^ I dKe 1Kr " K (4l) 

0 J — °0 

and using the formula for the Fourier transform of a normal frequency 
function and Weber's first exponential formula for the Bessel function 
integral, gives 


E(x,k) 




1/2 


exp 


l Q^ -f r c 
2t] 


(42) 


a normal distribution in both frequency and direction with rms values 
. Thus, the rms deviation in angle varies with the Mach number of 
the fluctuations and with the square root of both the mean eddy size 
and. the distance z. In addition to these, the relative frequency 
deviation also depends on the square root of the eddy frequency but is 
Independent of the eddy size. 

If the wave number part of the probability distribution equation (42) 
is now interpreted as the power spectrum of a wave, then it is possible 
to calculate the autocorrelation function of the signal at any position 
z. The normalized autocorrelation function is then 

' ,0 ° 

E(x,k)dk (43) 

0 

and for the above spectrum equation (42), 

\1/(t) = cos 03 o t exp ^-b^(a 0 T)^j (44) 

where b 2 = v^T) = M 2 n^Ar(9/2) z /l 92 j 2 c, and tie assumed relative "spread 11 
of the spectrum is slight. At z = 0, there is the cosine correlation 
function of a pure sine wave. For distances z, there is a loss in 
correlation for large values of delay t. 


?(t) 


!(x,k)cos u)t dk. 




CONCLUDING REMARKS 


It has been shown how the nature of a plane monochromatic sound 
wave may be changed by passage through a turbulent atmosphere. Perhaps 
it is well to pause and review the assumptions and then see how useful 
the results may be. 

The Lighthill and Kraichnan scattering formulae are valid when the 
Mach number of the turbulence is low (« l) . In addition, it will gen- 
erally be true that the sound particle velocity will be small compared 
with the turbulent fluctuation velocities. This means that the sound- 
turbulence interaction is assumed larger than the nonlinear acoustic 
effects. The sound generated by atmospheric turbulence usually will 
have frequencies much lower than those which are of interest in propa- 
gating through the atmosphere. The interaction or scattering phenomenon 
may then be considered to predominate in the audible frequency range. 

In evaluating the integrals in the appendices (C, D, E, and F) the 
assumption has been made that k 0 L » 1 , which means that the acoustic 
wavelength is small compared with the turbulent eddy size. This is not 
too severe a restriction since the scale of atmospheric turbulence is 
usually fairly large, of the order of the height above ground. In addi- 
tion, however, it was required that the scattering volume have dimensions 
large compared to L. This means that propagation distances z should 
be large compared to L, the scale of turbulence. 

With these restrictions, then, the transport-diffusion treatment of 
propagation In an isotropic turbulent atmosphere can predict the frequency 
spectrum and angular distribution of the sound wave fronts in a probabi- 
listic sense. The frequency distribution is interpreted as a power spec- 
trum and the autocorrelation of the wave form is obtained by Fourier 
inversion. Unfortunately the phase information is lost in the analysis. 

It is probable that the most significant failing of this work is 
its omission of the effects of thermal scattering which must occur when 
atmospheric turbulence mixes regions of different temperature. Unfor- 
tunately, the measurement of rapid temperature fluctuations in the atmos- 
phere is still in a rather rudimentary state. 


Massachusetts Institute of Technology, 
Cambridge, Mass., March 26, 1958. 



l6 


APPENDIX A 

RELATION BETWEEN SPECTRUM AND CROSS SECTION 


In reference 8, the power spectrum of the acoustic condensation £ 
scattered into the direction n’ from h is given by 

P(n,ca) = 1 5c | 2 |x(m)| 2 = bn ~ (S • S') 2 U n (K,D) 2 (Al) W 

2 

when K = k’ - k, Q = m' - m, and | x| is the distance from the scat- 7 

tering region. The symmetrical time transform of | is X(x,m) . The 
energy contained in a frequency interval dm' and a cone of solid angle 
d<J>' is thus 

P 

i P(n',m')p c 5 d$'dm' = ± P(S ' ,m’ dk ' (A2) 

Since by definition of the differential scattering cross section this 
must also be I 0 cr(k/k ' ) dk ' , it must follow ttat 

o(k/k ' ) = 2n ®^'(n • n' ) 2 U n (K,Sl)| 2 (A3) 



3T 
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APPENDIX B 
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SPECTRUM FOR ISOTROPIC TURBULENCE 


From equation (8), 




1 

(2k ) 4 


dx dx'dt 


iK(x-x ' )+ift(t-t ' 

CL o 6 



(Bl) 

where 


a = x - x 1 


t = t - t 1 


and 


Now, setting 




av 


® iJ (K,n) h 


u i u j 


av 


VT ( 2it ) ^ 

then it must follow by inversion that 


f io dTe- l(K - 3 - flT )R 1J (5,T) 


(B2) 


%j(5,0) = JdKe xK ' q J da3. i(3 (K,Q) 


(S3) 


where the quantity in brackets is the <J> i(j (K) introduced by Batchelor 
( re l* 9 , equation 2.4.3) • Accordingly, using reference 9 > 


/ wj 

= ® ± j (K) 

-00 

= - % k j ) 


8nK 


:(Aij - K iKj) r w(K,n)an (b4) 
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for Isotropic turbulence. The association may then be made 





(B5) 


whence follows 



(B6) 


where 0 is the angle between k and k T . 
unit volume and per second, set V = T = 1 


To seek the scattering per 
in equation (B6) . 


W 

1 

2 

7 
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APPENDIX C 


COSINE DRIFT COEFFICIENT 


Select the direction of k as a polar axis with an azimuthal 
angle y and a polar angle 0 with respect to the z-axis. The polar 
and azimuthal angles of k' with respect to k are 0 and |3, respec- 
tively. Using the identity 
W 

1 cos 6' = cos 6 cos 0 + sin 0 sin 9 cos(p - ty) (Cl) 

2 

7 the integral can he written as 

(y - l) - ^ jf 1 dxx5(l - 2x 2) 2 (1 - x 2)e- 8akx y' 2 ’ t d P {co S 6 (-2x 2 ) + 

sin 9 cos (P - (C2) 

where x = sin 0/2. The reader should consult Lighthill (ref. 7) for a 
discussion of the integral over 0 or x. For typical spatial turbu- 
lence spectra, assuming k/k Q » 1, the exponential diminishes very rapidly, 
and the upper limit of integration may be allowed to extend to infinity. 

The integration over p becomes -^iTtx^cos 0, and the integration over x 
is then 

/l • - s)„ « - | fj dxxV^ 1 " 

or 

r S m2 <«> 
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APPENDIX D 
COSINE VARIANCE 


write 


Using expression (Cl) and the spectrum of equation ( 15 ), one 


can 


(U' - l) 2 ) a = f Q dxx5 ( 1 - 2x 2 ) 2 (l - x 2 ) e J dp f_ 2x cos 0 

2 x(l - x 2 ) ' sin 0 cos(p - <|r) 2 J (d) 

The integral over (3 contains 8 itx 1 | cos 2 e and 2«x 2 (l - x 2 )sin 2 0. 

The x^ terms may be ignored compared to the x 2 terms (since the 
interest is in the range x « l) and this yields 


!*■ - |)2 X 


8kiAl0 


(l - ? 2 ) r x5(l - 2 x 2 ) 2 ( 1 - x 2 )e- 8akx (D2) 
J 0 


or 


ti * -»*)„= 1 2 ) * 

\ /o 52 a 


(D3) 


w 

1 

2 

7 
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APPENDIX E 

FREQUENCY DRIFT 


to 


From equations (26) and (22), the expression for 



is pro 


W 

1 

2 

7 


1 r °°( v 2 + p 2j -l^-v(p 2 -bpx 2 +kx 2 ) 1//2 


P dp 


where x = sin($/2) and A' = A/uj. The branch points of the intef 

are located at 


pi 

2 


hx^ ± \ l6x^ - l6x 2 + t 


0 ± 10/2 

±12 sin £ e y/ ~ ±12 sin (0/2) 


since the scattering angle is small. This makes ( 'p ^ proper tiona 

.1 -p(p 2 +4x 2 ) X/ ^ 2 


+ P 2 ) e 


/ (v 

J -00 

Hence /n\ ~ 0 

\ /o 


P dp 


■■parti 

(El) 

’"rand 

(E2) 

,1 to 

(E5) 


which must vanish. 
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APPENDIX F 


FREQUENCY VARIANCE 


Again, from equations (22) and ( 29 ), ^p 2 ^ Can ex P resse< ^ as 


8napAk 


7 r 1 


0 


dxx3(l - x 2 )(l - 2x 2 )V fakx T dp( 


x 2 2 2 
A + cu p 


-1 - ^p 2 

x x M 


(Fl) 


To evaluate this, calculate 

, 2„2 


f ape-^V + a , 2 p 2 )' 

-00 


where a 2 - ak/x. The integrand diminishes with p under the primary 
influence of the denominator. Thus the exponential is treated as a 
perturbation on this and written as follows 


2 2 


1 - a 2 p 2 


1 4- a 2 P 2 


assuming convergence for large p. Thus it is that one integral takes 
the form 



it 

A ( to + a A ) 


Now, 


/ 


p 2 dp 


2 2 

,-a p 


A 2 + a> 2 p 2 


da 2 A(a> + aA) 


RS « 

2auf 


(F2) 


if it is assumed that 
integration over x. 


03 


— » a. Thus, putting this result into the 
A 

the result becomes 


W 

1 

2 

7 



23 



4« 2 ap7qJ f 1 xT/^e" 80 -^ 

/ . xl/2 2 J 0 
(ak) ' a) c u 


= \f"2n 2 r(9/2) A m 2 
192 c 


(F3) 


W 

1 

2 

7 
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